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We study stationary and dynamical properties of the many-body Landau-Zener dynamics of a Bose 
quantum fluid confined in two coupled one-dimensional chains, using a many-body generalization 
recently reported [Y.-A. Chen et al.], within the decoupling approximation and the one-level band 
scheme. The energy spectrum evidences the structure of the avoided level crossings as a function 
of the on-site inter particle interaction strength. On the dynamical side, a phase diagram of the 
transfer efficiency across ground- state and inverse sweeps is presented. A totally different scenario 
with respect to the original single-particle Landau-Zener scheme is found for ground-state sweeps, in 
which a breakdown of the adiabatic region emerges as the sweep rate decreases. On the contrary, the 
transfer efficiency across inverse sweeps reveals consistent results with the single-particle Landau- 
Zener predictions. In the strong coupling regime, we find that there is a critical value of the on-site 
interaction for which the transfer of particles starts to vanish independently of the sweep rate. Our 
results are in qualitative agreement with those of the experimental counterpart. 



I. INTRODUCTION. 

The original Landau-Zener (LZ) problem involves a 
two-level system whose energy separation varies as a lin- 
ear function of time. These two states can be associ- 
ated to the lowest states of a single-particle in a double- 
well potential. As it is well known [H [2], the transi- 
tion probability between the two energy states is the 
well known LZ formula Plz = exp(— 27rJ^/to) written 
in terms of the tunnelling-coupling parameter J and the 
time-dependent detunning A = at, accounting for the 
probability that a particle initially in the bottom well at 
t = —oo reaches at t = oo the opposite well. Recent 
experiments have demonstrated the many body gener- 
alization of the LZ phenomena by loading an ultracold 
Bose gas into a pair of coupled one-dimensional chains 
where a controllable inter-chain sweep leads to the ob- 
servation of both, avoided crossings and breakdown of 
adiabatic inter-chain transfer [3 . 

Before the experimental realization of the many body 
LZ generalization, several theoretical descriptions of a 
Bose-Einstein condensate (BEC) confined in a double- 
well potential have addressed the LZ problem within two 
intrinsically different approaches. We make reference to 
the two-mode mean-field and Bose-Hubbard schemes in- 
herited from the Gross-Pitaevskii and full quantum ap- 
proaches [IHS]. As a result of incorporating the linear 
variation in time between the two levels, all of those 
treatments suggested a breakdown of the adiabatic limit, 
that is, that the LZ transition probability does not vanish 
even in the adiabatic limit. For example, the non-linear 
two-level system shows that the mean-field interactions 
among particles tend to increase the tunnelling proba- 
bility and that there exists a critical value of the inter- 
action strength beyond which the transition probability 
becomes nonzero even in the adiabatic limit f5^. On the 



other hand, a stationary phase approximation leads to 
a characteristic scaling or power law for the critical be- 
havior that occurs as the nonlinear parameter equals the 
gap of avoided crossing energy levels [8 . Regarding the 
asymmetric LZ tunnelling in a periodic potential, Jona- 
Lasinio et. al. [9] found that the tunnelling rates for the 
two directions of tunnelling are not the same. Tunnelling 
from the ground state to the excited state is enhanced by 
the nonlinearity, where as in the opposite direction it is 
suppressed. Even more, a LZ formula has been derived 
for the two mode many-particle scenario [7 . 

Here, we are devoted to the description of the many 
body LZ generalization described in the first paragraph. 
That is, we are interested in describing a BEC confined in 
two-coupled ID chains where the potential depths, defin- 
ing the lattice sites, are linearly modified in time such 
that a sweep from an initial energy difference —A to a 
final one A is achieved. To analyze such a system, we 
shall use a model Hamiltonian that incorporates both, 
the full quantum frame of the Bose-Hubbard Hamilto- 
nian [10 and a decoupling approximation scheme [TTJ[T2] 
that involves certain superfiuid order parameters. These 
are i/j^ = {h\^^) = {bi^u), where bl^^ and bi^^, v G {L,K\ 
create and annihilate particles in lattice sites i and side 
chain v G L^R. The average is the expectation value in 
the investigated states. A previous theoretical work has 
employed the time-dependent density matrix renormal- 
ization group method [13 where, with a small number 
of links in the chains, qualitative agreement with the ex- 
periment was found. Our model simulates two coupled 
infinite chains with a local site Hamiltonian written in 
terms of the on-site non-perturbative inter-particle inter- 
action strength [/, the time-dependent detunning A and 
the intra- and inter-chain coupling energies j" and J^. 
The model satisfactorily describes qualitatively the main 
features of the experiment: the transfer efficiency across 



ground- state sweeps, that in turn exhibits a breakdown of 
adiabaticity as the sweep rate decreases, and the transfer 
efficiency across inverse sweeps, that actually show con- 
sistent results with the single-particle Landau-Zener sce- 
nario. Regarding the single-particle LZ result the model 
Hamiltonian allows us to recover the result for a single 
dimer. In addition, since the number of particles per lat- 
tice sites in the actual experiments is not too large, it is 
possible to perform accurate numerical simulations of the 
system. Such a analysis gives rise to a richer spectrum 
with respect to the single dimer considered in [3]. 

The manuscript is organized in VI sections. In sec- 
tion II we introduce the model Hamiltonian that pro- 
vides a suitable description to the experiment. We also 
justify the use of the decoupling approximation and the 
one band scheme within the full quantum and mean- 
field approximation schemes. In section III we charac- 
terize the stationary states as a function of the both, the 
inter-particle interaction strength U and the the param- 
eter characterizing the tilt A. In particular, we concen- 
trate in determining the energy spectrum for the model 
Hamiltonian and the behaviour of the superfluid order 
parameters tp. An exhaustive exploration of the param- 
eter space allows us to encode the dynamical behavior 
of the of ground-state and inverse sweeps in a phase di- 
agram. The latter are presented in sections IV and V 
respectively. Finally, in section VI a summary of the 
main results is given. 



II. THE MODEL, DECOUPLING 
APPROXIMATION IN THE ONE-LEVEL BAND 
PICTURE 

The physical system that we want to describe is a fi- 
nite collection of interacting Bose atoms initially placed 
in one of the sides of a composed pairwise coupled lattice, 
see Fig. [ij and its subsequent evolution when the pair- 
wise potential depths are linearly modified in time. As 
schematically shown, the transport of particles in x and 
z directions occurs as a consequence of the intrinsic tun- 
nelling among the pairwise coupled sites, but it is also 
affected by the external variation of the well potential 
depths with time in the form of a sweep from an initial 
energy difference = — A to a final one Aj = A. The 
parameters responsible for such a transport are the on- 
site inter-particle interaction strength /7, the intra-chain 
coupling energy j" {z direction), the inter-chain coupling 
energy {x direction) and the time-dependent param- 
eter characterizing the tilt between the two wells of each 
site in the chain, namely, the detunning parameter A. 

The purpose of the present paper is to fully charac- 
terize stationary and dynamical properties of the system 
described above. To proceed, we will consider the ex- 
perimental realization reported in Ref. ^ where the 
transfer of particles in x and z directions is studied in 
what it is called the ground-state and inverse sweeps, the 
former accounting for a sweep starting from the ground 
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FIG. 1. (color online) Schematic representation of the double- 
well trap chain potential. The coupling energies are intra- 
chain j" (z direction) and inter-chain (x direction). The 
time-dependent parameter A characterizes the tilt between 
the two wells of each site in the chain. 

state, namely, where all the particles initiate in the bot- 
tom well with lower energy, and the latter where the ini- 
tially filled sites are the ones with higher energy. This 
scheme leads us to answer several questions, on one side 
those related to the stationary properties, such as, the 
energy spectrum, and on the other side, aspects related 
to the dynamical behavior of the atomic population in 
the wells. Of special interest in this regard is the deter- 
mination of the phase diagram encoding the dependence 
of the transfer efficiency on the parameters characterizing 
the dynamics. 

Let us now turn our attention to the model Hamilto- 
nian of the system. Since the particles are not allowed 
to tunnel between left and right wells of different sites 
i, movements in x and z directions must be considered 
separately. Therefore, we consider the dynamics along 
the sites and their counterpart between the left and right 
chains. 

Regarding the two-well geometry and its intrinsic 
asymmetry supplied by the time-dependent parameter 
A, we should first argue about the level picture scheme 
where our system is appropriately described. Due to the 
variation of such a parameter not just the first band but 
excited levels might play a significant role. As estab- 
lished by Dounas-Frazer et. al. [10] the transport phe- 
nomena of a Bose-Einstein condensate in an asymmetric 
double well leads to three new energy scales with respect 
to the symmetric trap and one-band assumptions. These 
are the hopping, tilt and the energy gap between the 
first and second bands. The task of including at least a 
two-band picture (that is, four levels per site) leads to 
a Hilbert space of size {N + l)(Ar + 2){N + 3)/6 with 
N the number of atoms per site. However, since the LZ 
phenomena involves the presence of just two states, we 
should restrict ourselves to the first band of the time- 
varying tilted double-well potential. This can be done 
by appropriately taking the formal bounds where a one- 
band picture is well justified, which primarily assumes 
that the interaction energy is much smaller that the en- 
ergy band difference. In the appendix we present a study 



of the single-particle energy levels in a double- well poten- 
tial as a function of the tilt A. To fulfill the condition 
that the inter-particle energy interaction is much smaller 
than the energy gap between the first and second bands 
we shall confine ourselves to values of A for which the 
bands are well separated and thus the condition with the 
inter-particle interaction will be considered accordingly. 

Starting from the Bose-Hubbard Hamiltonian of an op- 
tical lattice [3 , one can write down the Hamiltonian of 
atoms in the two- well chain as, 

i i 
u=R,L <i,j> 

i 

= ^ E blAAA,'^ (1) 

u=R,L 

where the labels || and ± denote movements in z and 
X directions respectively. represents the tunnelling 
along the chains, while Tif- takes into account the local 
two-level structure at sites i. Hf corresponds to the on- 
site interaction energy. 

To consider the participation of the lattice sites, that is 
the transport in the z direction, we use the decoupling ap- 
proximation where order parameters ipi, = {h\^) = (bi^u), 
V G {L^R} are introduced to account for the expecta- 
tion value of the tunnelling of particles between left and 
right wells, disregarding the position of the site i within 
the chain. We consider negligible the average fiuctua- 
tions of creation and annihilation particle operators in 
sites i and j, such that, A6j^ = — ijjy ^ and 
A6^-,^ = bj.v-'^v ^ 0. This allows us to write [in [T2l[T4] . 
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FIG. 2. (color on-line) Schematic representation of the lattice 
sites in the chain of double- wells. 

and — and sides R and L. The Fock space elements 
can be written as ^i^^) where = 0, ...,A^, 

ni = 0,...,Ar-n*^ and n^^ = N-ni^-ni, being AT the 
total number of particles at each site i. Thus, our model 
Hamiltonian, representing the local site energy of an in- 
finite chain, corresponds effectively to N particles with 
three "internal" states. Consequently the Fock space at 
each lattice site scales with the number of particles as 
n = (A/' + l)(7V + 2)/2. 

It is important to notice that while the transport along 
z direction is accounted for in a mean-field fashion, the 
tunnelling among left and right sides {x direction) is fully 
quantum described. The assumption of the mean-field 
approach in this case is well justified in the sense that 
in the experimental situation near the centre of the trap, 
the average the number of particles per site on left and 
right sides, is the same along the chains. As we shall see 
in the next sections, this treatment allows us to study up 
to = 8 particles per site. 



III. STATIONARY STATES 



Thus, by substituting this expression into we arrive 
to the model Hamiltonian at each lattice site given by, 

^loc^.^ii ^ (^^(s,,+sg-^3)+H,^+^f. 

u=R,L 

(3) 

This model Hamiltonian represents the local site en- 
ergy of an infinite chain, where the interaction is non- 
pert urbative. We should emphasize that the fact of 
considering a mean-field approach to separate the ef- 
fective contributions by site leads us to the task of 
self-consistently determining the order parameters ^/^^y, 
u G {L^R} to specify the Hamiltonian ([3|. 

Regarding the Hilbert space size, one can see that the 
Fock space of the coupled chains can be described in 
terms of the occupation number of sites i and In Fig. 
[2] we illustrate the equivalence among both, sites i — (i-fl) 



In this section we study the stationary states of Hamil- 
tonian ^ for a given number of particles, as a function 
of the interaction strength U and the parameter char- 
acterizing the tilt A. In particular, we solve the time- 
independent Schrodinger equation for the model Hamil- 
tonian to determine the energy spectrum as a function 
of the parameters U and A. It is important to point 
out that in our calculations both, the parameters U 
and A are referred to the tunnelling coupling parame- 
ter J = jll = J^. Here and henceforth we use U = U/J 
and A = A/J. 

To determine the entire energy spectrum of Hamilto- 
nian ([3| we have to find the value of the order param- 
eters tjjjy {v = L^R). We proceed as follows. For fixed 
values of JJ and A we obtain by means of a variational 
procedure the values of iIjj, {v = L^R) that mini- 
mize the ground state energy Eq(U ^ A, ipL^i^n)- In other 
words, we search for the lowest eigenvalue of the Hamil- 
tonian hJ^^(/7, A^ipL^ ipn) by minimizing with respect to 
t/jjy {v = L^R). Once those optimal values have been de- 
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FIG. 3. (color on-line) Density plot of the optimal values of 
the order parameters ipL and ipR as functions of the dimen- 
sionless parameters U and A for the ground state of Hamil- 
tonian ([3| for iV = 6. 



termined, we calculate the entire spectrum by means of 
standard diagonalization methods. Then, for the same 
value of U we repeat the same process for each value of 
A 

In Fig. |3]we show a density plot with the optimal val- 
ues of the order parameters ipL (left) and ipR (right), as 
functions of U and A. These calculations performed for 
= 6, allow us to appreciate well defined boundaries 
for which the optimal values of the order parameters are 
identically zero. These are Mott insulating (MI) state- 
like regions. That is, according to the identification given 
in section [nj tjjj^ = {h\^^) = {hi^v)^ v G the re- 

gion delimited for those values of U and A indicates that 
there is no tunnelling and consequently a Mott Insulator 
phase can be nucleated since the average fluctuations in 
the population per site are negligible. Conversely, a su- 
perfluid phase is associated when the order parameters 
are different from zero, the superfluid component in the 
system being maximal in the diffuse area for IJ < 10. 
As expected, the indistinguishable character of left and 
right wells is confirmed from the symmetry exhibited in 
and i\)R. It is important to point out that the insulat- 
ing regions that appear for ~ and large | A | are finite 
size effects similar to what was found in small Bose- Fermi 
systems [16 , which in the limit of large N these regions 
disappear. 

For illustration purposes, in Fig. [4| we show the en- 
tire energy spectrum of HamiltonianTsl) for A/" = 3 and 
different values of IJ . The values of and A have 
been rescaled with respect to IJ as A^ = KjlJ and 

^ ^^/^' ^ = numbers on the right 

refer to the number of particles in the right well, while 
the circles indicate how tunnelling resonances are crossed 
one after another. The left figures ((a), (b) and (c)) are 
obtained by neglecting the transport along the z axis (see 
Fig.l). That is, they correspond to the energy spectrum 
of a single double-well. The right figures of the panel 



contains the energy spectrum considering the transport 
in both, parallel and perpendicular directions. From 
the full panel of Fig. [4] one can observe that although 
avoided level crossings are exhibited in the whole inter- 
val [— A^, A^] for every considered value of the level 
crossing structure is more visible as U grows. As can be 
appreciated, the results of our model (|3| when the trans- 
port along the z axis is neglected, coincide with those 
obtained for Bose-Hubbard Hamiltonians of one dimer 
only O [12 . In other words, if the on-site interaction 
energy U between the particles dominates the tunnelling 
coupling J {U = U / J » 1) ^ the LZ transitions split into 
N avoided level crossings [see Figj4](c)], each correspond- 
ing to the transfer of one particle. A different situation 
is observed when particles can move through x and z 
axis. The LZ transition for U = U/J » 1 split into 
2N — 1 avoided level crossings. Regarding the number 
of avoided crossing and its dependence with the number 
of particles N in the limit of non-interacting bosons, the 
single-particle LZ result is recovered for weaker interac- 
tions when the transport in x and z axis is considered. 
It is important to emphasize that the fact of considering 
the combined transport in the parallel and perpendicu- 
lar directions, in general gives rise to a richness spectrum 
with respect to the single dimer. Even more, leads to pre- 
dict additional tunnelling resonances with respect to the 
single dimer. 



IV. PHASE DIAGRAM 

The LZ dynamics from A = —aTmax to A = aTmax 
for a given value of a are obtained by evolving in time 
(A = at) the initial state |^(0)) associated to the ground 
state energy Eo{U^A). We have that, for a given value 
of a one should follow the evolution of the state where 
all the particles initiate in the left chain with A < 0, 
such that the filled chain is the one with lower energy. 
Experimentally, this process is identified as a ground- 
state sweep. The evolution is given by [17], 



(4) 



where H is the local Hamiltonian ([3| and h = 1 is 
assumed. The initial time to corresponds to Tmax = 
— Ao/a while the temporal step St that we select, St = 
0.2, produces qualitatively similar results than St < 0.2. 
It is important to emphasize th at t he minimization pro- 
cedure described above (section [TTT| ) is performed at each 
temporal step St to take into account the variation of 
A at each time interval St. Consequently, the ground 
state l^(ti)) is updated at each temporal step and serves 
as a seed for the subsequent time. The full dynamics, 
for a given rate a and inter particle interaction strength 
/7, ends when the time t^+i reaches its maximum value 
Tmax = Ao/a. To illustrate the evolution in time of the 
population on the right side, in Fig. [5|a,b,c) we plot 
nR{t) for various rates a, U = 0.75 and = 6. The 



5 




-4 -2 2 4 -4 -2 2 4 



FIG. 4. (color on-line) Energy spectrum = 
En/JU of Hamiltonian ([s]), for = 3 and ?7 = 

1 (a,d),10 (b,e),100 (c,f). The values of A in the x axis 
have been rescaled with respect to the dimensionless param- 
eter that is, = A/U. For [/ = 100 the numbers on the 
right refer to the number of particles in the right well, while 
the circles indicate the tunnelling resonances. 



i = t/Trr 



with 



time is rescaled with respect to 
Ao = 20. 

For A/" = 6 we investigated the many-body dynam- 
ics for sweep rates in the interval 27r/a G [0.1,10] and 
U G [0, 10] [19]. To condense the information of the LZ 
dynamics, as a function of the interaction strength U 
and the sweep rate a, we show in a phase diagram the 
final state that the system reaches as the initial state 
is evolved in the interval [-Tmax^Tmax]- In Fig. [5|d) 
we show a density plot with the normalized transfer effi- 
ciency on the right side at t = Tmax- To complementary 
visualize the dependence of ur on U in Fig. [sje) we 
present several contour plots of the normalized transfer 
efficiency as a function of the sweep rate a, in the interval 

il e [0,3]. 
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FIG. 5. (color on line) Instantaneous right evolu- 
tion riRit) from ground-state for various rates 2tt / a — 
0.1(a), 1.0(b), 2.0(c), U = 0.75 and AT = 6. The time is 
rescaled with respect to Tmax, t = t/Tmax, with Aq = 20. 
(c) Density plot of the normalized transfer efficiency ur at 
t = Tmax- (d) Contour plots of the normalized transfer effi- 
ciency for several values of the the dimensionless parameter 
U, see main text. 



From Fig. [5|d) one can observe how the influence of 
both, the parameter characterizing the interaction among 
particles and the sweep rate a, lead the system to well 
defined quantum phases from minimum to maximal val- 
ues of the transfer efficiency nn. As can be observed from 
this figure, the effect of having N > 1 induces the oppo- 
site behavior with respect to the original LZ phenomena 
when the sweep rate a decreases. That is, the transfer of 
particles does not vanish as the sweep rate tends to zero, 
on the contrary, when the sweep rate a takes smaller 
values, the total population initially placed on the left 
side is transferred almost completely to the right chain. 
In other words, the main effect of having an interacting 
many-body system is the breakdown of adiabaticity. 

Regarding the role of the on-site interactions, one can 
observe that the transfer efficiency diminishes as U in- 
creases. This result can be understood in the light of the 
well known prediction for bosonic Josephson junctions 
[20] where the transport of particles diminishes as the in- 
teraction strength increases. We observe that for a given 
value of the on-site inter-particle interaction the effect of 
decreasing the sweep rate a, is turned from detrimental 
into favourable for the transfer efficiency in the A/"— body 
LZ scheme \see Fig. [5](d) and (e) ]. 

To illustrate the dependence of the many-body dynam- 
ics on the initial value of the tilt — Aq and on the number 
of particles, we performed calculations for Aq = 10 and 
100 and N = 2 and 8 particles per site. The results 
are shown in the panel of Fig. |6] As above, we plot 
the normalized transfer efficiency in a density plot for 
27r/a G [0.1,10.0] and UN/Aq G [0,3]. The phase dia- 
grams in Fig. [ojevidence the impact of the initial tilt — Aq 
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FIG. 6. (color on-line) Density plot of the normalized transfer 
efficiency ur of ground-state sweeps as a function of dimen- 
sionless parameters U and a at t = Tmax = Ao/a for N — 2 
(a,b) and, N = 8 (c,d). — Aq is the initial value of the tilt. 
Parameters are: Ao = 10 (a,c) and Ao = 100 (b,d). 



on the final state that the system reaches after a com- 
plete sweep. The region of maximum transfer is reduced 
for either N = 2 and A/" = 8 as Aq goes from 10 to 100. 
That is, the region exhibiting breakdown of adiabatic- 
ity becomes larger as the value of time dependent energy 
bias grows. From our simulations we find that there ex- 
its a critical value of the on-site interaction Uc ~ Aq/A/" 
for which the atoms remain evenly distributed. Due to 
the relation with the order parameters [Fig. [s], we find 
that the transfer efficiency is very sensitive with respect 
to the variation of ipL and iJjr. As the system evolves 
in time the presence and the crossing between superfluid 
and MI phases, give rise to reach a saturation value for 
the of the atoms transfer to the right side of the chain 
(max(n,^)<0.5for /7>/7c). 



V. MANY-BODY DYNAMICS: 
GROUND-STATE VS INVERSE SWEEPS 

Following the experimental realization of LZ many- 
body dynamics we extend our study to investigate 
the dynamics across inverse sweeps. Namely, the initial 
condition corresponds to A > with the total particle 
population on the left side. The time evolution is given 
as before, by Q, and the same minimization procedure 
to determine the variational order parameters iJjl and 
ipR at each temporal step St must be carried out from 
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FIG. 7. (color on-line) Density plot of normalized transfer 
efficiency ur of inverse sweeps as a function of dimensionless 
parameters U and a for 27r/a G [0.1,10.0] and UN/ Aq G 
[0,3]. Parameters are: Aq = 10 (a), Aq = 100 (b), = 2. 



A = aTrr 



Ao to A 



-aTrr 



-Aq for a given 



sweep rate a. 

To obtain the phase diagram for the inverse sweep 
scheme we explore sweep rates in the interval 27r/a G 
[0.1, 10.0] for N = 2, and set Aq = 10 and Aq = 100. For 
comparison purposes with the results presented in Fig. 
|6) we plot in Fig. [7| the phase diagram of ur in terms of 
UN/ Aq. In this case consistent results with the single- 
particle LZ scenario are found: the maximum of transfer 
efficiency is observed as the sweep rate increases. The 
effect of increasing the value of the tilt Aq from 10 to 
100 gives rise to a reduction of the region of maximum 
transfer. That is, it leads to similar qualitative features 
than those found for the ground state sweeps. The role 
of the inter particle interaction strength is perceived as 
in the ground-state scheme, that is, increasing the value 
of U destroys localisation between left and right chains 
{riR ^riL < 0.5). 

To compare the dynamical behavior between the 
ground-state and inverse sweeps we study, as in the ex- 
perimental counterpart, the transferred population as a 
function of the final value of the tilt Aj, for N = 2, a 
fixed sweep rate a and fixed interaction U. As before, 
the initial state for the ground-state (inverse) sweep cor- 
responds to the total particle population localized in the 
bottom well with lower (higher) energy. However, it is 
important to point out that in this case the sweeps are 
carried out from A^; to Aj, that is, the sweeps are not 

As ini- 
20 and 



IV 



accomplished from — |A^| to |A^| as in section 
tial condition for each sweep we consider A^ 
A^ = —20 for ground-state and inverse sweeps respec- 
tively. Fig. [s] corresponds to ground-state (a,c) [ inverse 
(b,d)] sweeps. In Fig. [s] (a) [(b)] we plot ur for ground- 
state [inverse] sweeps as a function of Af for U = 1.5 
and several values of 2ix /ol. In Fig. [8](c) [(d)] we plot ur 
for ground-state [inverse] sweeps as a function of Aj for 
27r/a = 2.0 [1.0] and U G [0, 10]. As can be observed from 
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FIG. 8. (color on-line) Normalized transfer efficiency tir for 
ground- state (a) [inverse (b)] sweeps as a function of the final 
tilt A/, for [/ = 1.5 and 2ti / a G [0.1, 3.0]. Normalized transfer 
efficiency ur as a function of the final tilt A/, for 2ti / a — 2.0, 
[2n/a = 1.0] in the ground-state (c) [inverse (d)] case and 
tf G [0, 10]. The number of particles is N — 2. Dimensionless 
parameters are: ur^ A/, U and a. 

these figures the ground state sweeps show that the nor- 
malized population transferred increases from zero to a 
maximum max(ni^) < 1 for 27r/a < 1 while the behavior 
for the inverse sweeps shows that the population trans- 
ferred is partial, reaching a maximum value below one. 
The effect of increasing interaction for a fixed sweep rate, 
is to suppress the maximum population transfer. Com- 
parison with the experimental results show qualitative 
agreement. See Fig. 3 of Ref. [3 . 



VI. FINAL REMARKS 

We have studied the dynamical and stationary proper- 
ties of an interacting Bose quantum fluid confined in two 
coupled one-dimensional chains. Such a system has been 
realized experimentally [3 and represents a generaliza- 
tion of the single-particle Landau-Zener dynamics. The 
competition terms that determine the evolution in time 
from a given initial state depend on the inter particle 
interaction strength [/, the intra-chain coupling energy 
jll {z direction), the inter-chain coupling energy {x 
direction) and the time-dependent parameter character- 
izing the tilt between the two wells of each site in the 
chain, that is, the detunning parameter A. 

To investigate the dynamics of the generalized LZ real- 
ization, we work in the decoupling approximation and the 
one- level band scheme, and propose a model lattice site 



Hamiltonian that represents two coupled infinite chains 
(|3|. Our full many-body model is written in terms of 
the order parameters ifj^^ v ^ {L^R\ that account for 
the superfluid component of the particle population. We 
exhaustively explore the space of parameters {/7, A} to 
determine by variational means the order parameters ipy^ 
V G {L, i?} for the ground state energy Eo{U^ A^^jj^^ i^n)- 
Then, we diagonalized the Hamiltonian to obtain the en- 
ergy spectrum. Concerning the stationary state prop- 
erties, since the number of particles per lattice sites in 
the actual experiments is not too large, we were able to 
perform accurate numerical simulations of the system to 
fully characterize the entire energy spectrum. Such cal- 
culations lead to predict a richer spectrum with respect 
to that performed for a single dimer, that is, when the 
movement along z axis is neglected [3 . We found that 
as U grows in magnitude, a more complex structure of 
the energy spectrum with respect to the double- well sys- 
tem is developed in correlation with the suppression of 
the superfluid component in the system, thus promoting 
self-trapping behavior. In addition our model allow us to 
recover the single-particle LZ result for a single dimer. 

Regarding the dynamical evolution of the LZ many 
body generalization, we concentrate in studying ground- 
state and inverse sweeps in order to compare the pre- 
dictions of our model with the experimental realization. 
The ground-state and inverse sweeps dynamics account 
for a sweep from a state where all the particles initi- 
ate in the bottom well with lower energy (left well) and 
the counterpart where the filled sites are the ones with 
higher energy, respectively. In a phase diagram in terms 
of the interaction U and the sweep rate a we condense 
the normalized final transfer of atoms. A totally differ- 
ent scenario with respect to the original single-particle 
Landau-Zener scheme was found for ground-state sweeps. 
That is, the transfer of particles does not vanish as the 
sweep rate tends to zero, on the contrary, when the sweep 
rate a takes smaller values, the total population initially 
placed on the left side is transferred almost completely to 
the right chain. A different situation occurs in the study 
of the transfer efficiency for inverse sweeps, this reveals 
consistent results with the single-particle Landau-Zener 
phenomena predictions. As the system evolves in time 
depending on the interaction, the system can transit dif- 
ferent phases (superfluid or insulating) multiple times, 
see Fig. [3l We identified a a critical value of the interac- 
tion Uc, for which U > Uc, the system is in a crossover 
from insulating states {ur ^ ul) with iIjr^l = to a 
superfluid states {nji ^ ul ^ 0.5) with i/jr^l 7^ . It 
is important to emphasize that all of our results are in 
qualitative agreement with the experimental realization. 
In addition, our analysis provides insight to the region 
where the interaction is non-perturvative. 

This work was partially supported by grants PAPIT 
IN108812-2 and IN116110 DGAPA (UNAM). Compu- 
tational facilities by DGSCA-UNAM are also acknowl- 
edged. We would like to thank V. Gal for useful dis- 
cussions regarding the implementation of computation 
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FIG. 9. (color on-line) Energy spectrum of a single-particle in 
a double well potential. The asym metry is determined by the 
value of the parameter c in (A.l). In these plots c = 0,0.04 
and 0.08. 



algorithms. 



Appendix: Single particle in a tilted double well 

To establish appropriate bounds for the values of the 
parameters that determine the dynamics of the Bose gas 
in the coupled one-dimensional chains, that is, their prac- 
tical use in relation to the proposed model Hamiltonian 
(|3|, we analyze here the eigensystem of a single particle 
confined in an asymmetric double well. Although a pre- 
vious study has provided quantitative criteria regarding 
the energy level scheme where a many-particle system 
can be properly described for those potentials p^, the 
intrinsic nature of the Landau-Zener phenomena, as well 
as, the proposed model require to justify the plausibility 
of the chosen parameters. 

Let us consider the potential 



where d = 0.05 and the parameter c incorporates the 
asymmetry. In Fig. [9] we show the energy spectrum, as 
well as, the single-particle wave functions corresponding 
to different values of c. From this figure we observe sev- 
eral features. As expected, the energy levels are arranged 
in bands of two levels. The energy difference between 
bands becomes similar to the energy difference intra- 
bands as the asymmetry in the potential grows. That 
is, as the asymmetry in the double-well potential, repre- 
sented by the parameter A = at changes with time the 
energy spacing between bands starts to match the energy 
difference among levels in the bands. This fact imposes 
a restriction on the possible values of the A parameter. 
In principle, one could think that the suitable range for 
values of A is severely limited. However, we recall that 
in our calculations all the parameters are referred to the 
tunnelling coupling parameter J = j" = J^, which is 
related to the energy difference of the two first energy 
levels at t = (c = 0) as J = (ei — eo)/2 pT], being 
eo and ei the ground state and the first excited state 
energies of the single particle. Therefore, by restricting 
ourselves to those values ofA = A/J « hw/J and 
U = U / J « hjj / J ^ being hhj ^ (es + 62 — ei — eo)/2 the 
energy difference among bands, we can validate the one- 
band scheme to determine the dynamics of the proposed 
mod el Ha miltonian. For the chosen double-well poten- 
tial (A.l) we found hw/J ~ 1.2 x 10^ and consequently 

^max/ J ^ 1-2 X 10^. Thus, for values of A < Arnax/ J 

the choice of parameters in the simulations and our model 
are well justified. 



REFERENCES 



V{x) 



(A.l) 



[1] L.D. Landau, Phys. Z. Sowjetunion 2, 46 (1932). 

[2] CM. Zener, Proc. R. Soc. London A 137, 696 (1932). 

[3] Y. Chen, S.D. Ruber, S. Trotsky, 1. Bloch and E. Altman, 

Nature Phys. 7, 61-67 (2011). 
[4] O. Zobay and B.M. Garraway, Phys. Rev. A 61, 033603 

(2000). 

[5] B. Wu and Q. Niu, Phys. Rev. A 61, 023402 (2000). 

[6] B. Wu and Q. Niu, New. J. Phys. 5, 104 (2003). 

[7] D. Witthaut, E.M. Graefe, H.J. Korsch, Phys. Rev. A 

73, 063609 (2006). 
[8] J. Liu, L. Fu, B. Y. Ou, S. G. Chen, D. I. Choi, B. Wu, 

Q. Niu, Phys. Rev. A 66, 023404 (2002). 
[9] M. Jona-Lasinio, O. Morsch, M. Cristiani, N. Malossi, J. 
H. Miiller, E. Courtade, M. Anderlini, and E. Arimondo, 
Phys. Rev. Lett. 91, 230406 (2003). 
[10] D. R. Dounas-Frazer, A. M. Hermundstad, and L. D. 

Carr, Phys. Rev. Lett. 99, 200402 (2007). 
[11] D. van Oosten, P. van der Straten, and H.T. C. Stoof, 

Phys. Rev. A 63, 053601 (2001). 
[12] M.P.A. Fisher, P.B. Weichman, G. Grinstein and D.S. 



Fisher, Phys. Rev. B 40, 546 (1989). 

[13] C. Kasztelan, S. Trotzky, Y.-A. Chen, 1. Bloch, 1. P. Mc- 
Culloch, U. Schollwock and G. Orso, Phys. Rev. Lett. 
106, 155302 (2011). 

[14] I. Danshita, J.E. Williams, C.A.R. Sa de Melo, and C.W. 
Clark, Phys. Rev. A 76, 043606 (2007). 

[15] We use the Nelder-Mead simplex method to find the min- 
imum of the surface defined by U/J and A/ J for the 
lowest eigenvalue of the Hamiltonian (|3| wth fixed N . 
Numerical Recipes in C. The Art of Scientific Comput- 
ing, 2nd Edition, 1992. 

[16] S. F. Caballero-Benitez and E. A. Ostrosvskaya, J. Phys. 
B: At. Mol. Opt. Phys. 44, 135301 (2011). 

[17] The approximation of the time evolution is accurate to 

[18] T. Venumadhav, Masudul Haque and R. Moessner, Phys. 

Rev. B 81, 054305 (2010). 
[19] The interval of study in 2Ti/a and U is 10/2^, that IS, we 

explore a grid of size 2^ x 2^. 
[20] S. F. Caballero Benitez, V. Romero-Rochm and R. Pare- 



9 



des, J. Phys. B: At. Mol. Opt. Phys. 43, 115301 (2010). 
[21] R. Paredes, Phys. Rev. A 73, 033616 (2006). 



